Violation of Porod law in a freely cooling granular gas in one dimension 
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We study a model of freely cooling inelastic granular gas in one dimension, with a restitution 
coefficient which approaches the elastic limit below a relative velocity scale 5. While at early times 
{t <g S~^) the gas behaves as a completely inelastic sticky gas conforming to predictions of earlier 
studies, at late times {t 3> 5~^) it exhibits a new fluctuation dominated phase ordering state. We 
find distinct scaling behavior for the (i) density distribution function, (ii) occupied and empty gap 
distribution functions, (iii) the density structure function and (iv) the velocity structure function, 
as compared to the completely inelastic sticky gas. The spatial structure functions (iii) and (iv) 
violate the Porod law. Within a mean-field approximation, the exponents describing the structure 
functions are related to those describing the spatial gap distribution functions. 

PACS numbers: 47.70.Nd, 64.75.+g, 05.70.Ln, 45.70.Mg 



Flowini 
nomena 



xanular media exhibit varied physical phe- 
3|. A simple, well studied model that cap- 
tures many features of flowing granular systems is agas 
of particles undergoing inelastic collisions [1, H, 0, 0, [a S 
The system may be externally driven or allowed 
to cool freely. However, real granular gases have elas- 
tic collisions when relative velocity of particles approach 
zero [llj. Thus, a reahstic model of cooHng granular gas 
should have a relative velocity dependent restitution co- 
efficient pjj. In this paper we focus on such a model. 

In general, a system freely relaxing to an ordered state 
has a macroscopic length scale C{t) increasing with time 



t IJI- In addition, for usual phase ordering systems, the 



presence of a dominant i2(t) results in a robust scaling 
law called the Porod law [13) ■ For scalar order pa- 
rameters, the Porod law states that the scaled structure 
function S/ ~ {kC)~^ for large kC, with 6* = 2 in one 
dimension. Contrary to this clean scenario, in certain 
driven systems [3, E, [3 12 1 1 an unusual phase order- 
ing was observed. These systems have a macroscopic 
coarsening length scale C but the domain lengths have a 
power-law distribution with a large negative power. This 
leads to a violation of the Porod law with 0^2. We shall 
refer to such systems as fluctuation-dominated phase or- 
dering (FDPO) systems. We demonstrate below that a 
coarsening granular gas too has such an unusual FDPO 
state. 

The issue of an inelastic gas showing coarsening and 
phase ordering has been addressed in many earlier pub- 
lications on the subject [E 0, S B K 111, Ell- It is 
now well-known that both the freely cooling inelastic gas 
(0 < r < 1) and the sticky gas (r = 0) undergo coarsen- 
ing with a growing length scale C{t) ~ t^/^ , with z — 3/2 
in one dimension S S|- The sticky gas problem in 
one dimension can be solved exactly and is known to be 
equivalent to the inviscid Burgers equation [1, . From 
this equivalence, the structure functions of the sticky gas 
can be inferred to obey the Porod law ji^l . 



Numerical studies have tried to relate the behavior of 
the inelastic gas, to the sticky gas. It was shown that at 
large times the decay of total energy of the inelastic gas is 
identical to the sticky gas Moreover, other quantities 
like C{t) and velocity distribution function have the same 
scaling form as that of the sticky gas. This suggested 
that for any deviation from the elastic limit, the large 
time scaling behavior crosses over to that of the sticky 
gas and that the underlying continuum equation is the 
inviscid Burgers equation [6|. 

In this paper, we show that if a granular gas is mod- 
eled as a gas of particles having a velocity dependent 
restitution coefficient, interesting new physics appear. In 
the late time regime, the equivalence with the sticky gas 
breaks down and the system exhibits a FDPO state. For 
our granular gas model, velocity dependence of the resti- 
tution coefficient is chosen to be: 



r = (1 - ro)exp {-\v,c\/5\") -f tq. 



(1) 



For relative velocity v-ccx <C (5, r — > 1, and for Urci 3> 5, 
r — !■ ro < 1 , mimicking the experimental scenario 
The parameter a determines how sharply the crossover 
from ro to 1 happens across the crossover scale 5. While 
experiments suggest a wide range of values for cr, 
kinetic theory studies estimate a — 1/b [21j. We note 
that taking first the limit cr — *■ cxd and then (5 — > 0, this 
model becomes the same as the model studied in Ref. 0] . 

Our main result is that the cooling granular gas has a 
time scale ti ~ , such that the density distribution 
function, and various spatial distribution functions show 
a complete change of behavior across it. Yet at the same 
time, the total energy E{t) decay as t~'^/'^ throughout, 
without any signature of change across the timescale ti. 
For t <C ti, we find that the granular gas behaves as a 
sticky gas (as in j6]). But, for late times ti <^ t <^ t2 
(where t2 ~ the phase ordering is distinct from 

the sticky gas. In particular, the density-density and 
velocity- velocity structure functions show violation of the 
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FIG. 1: The variation of (p^) with time t is shown for different 

values of S when a = 3.0. The dashed Une corresponds to the FIG. 2: The variation of Pi{p, t) with p is shown for different 
sticky gas. Inset: The curves collapse when t is scaled by . times. 



Porod law. We note that the scale ^2, beyond which all 
collisions are elastic and E{t) stops decreasing, is eas- 
ily understood in terms of the velocity scaling law Q - 

v{t2) ^ ^2 ^^"^ ^ ^ implying t2 ^ S~'^. On the contrary, 
the interesting scale ti that we find is much smaller. 

We now define our model more precisely. We consider 
N point particles of equal mass on a ring of length L. 
Initially, the particles are distributed randomly in space 
with their velocities drawn from a normal distribution. 
The particles undergo inelastic, momentum conserving 
collisions such that when two particles with velocity Ui 
and Uj collide, the final velocities u[ and u'^ are given by: 

We define coarse-grained densities and velocities for the 
granular gas as follows At any point of time the 
system is divided into N equally sized spatial boxes. The 
total number of particles in the ith {i — 1,2, ... , N) box 
defines the mass density pi. The velocity Vi is defined as 
the sum of the velocities of the particles in box i. For the 
sticky gas it suffices to talk about distributions of masses 
and velocities of individual particles. 

We have done an event driven molecular dynamics sim- 
ulation for system sizes L, ranging from 20000 — 50000 
(in units of inter particle spacing). The particle density 
is set to 1 throughout. The system was evolved up to 
time t — 32000 — 64000 (in units of initial mean collision 
time). For these times C{t) <C L. Simulations were done 
for different a (namely 3, 4, 5, 10 and oo), 5 (namely 
0.001 - 0.01), and ro (namely 0.1, 0.5 and 0.8) values. 
There was no qualitative difference found for the various 
sets of these parameter values. Hence, we choose a spe- 
cific set for the data presented below - — 0.5, a — oo, 
5 = 0.01 (unless otherwise mentioned). Whenever there 
is a quantitative dependence on a, we mention it. 



The existence of the time scale ti can be seen by exam- 
ining (p^), where p is the density. In Fig. [1] we show the 
time dependence of (p^) for systems with same cr = 3 but 
different 5 values. At early times {p^) ^ t^/"^ as for the 
sticky gas. The departure from the sticky gas curve hap- 
pens at a timescale ti which increases with decreasing 6. 
The curves collapse when t is scaled by 6^^ (see inset). 
We have checked that this dependence is independent of 
a, i.e., ti for all a. 

Currently we do not have a deeper understanding of 
the timescale ti. Intriguingly, we find no signature of ti 
in the decay of the total energy E{t) ^ t^^^"^, which is re- 
lated to the second moment of the velocity distribution. 
In this paper, we concentrate on the density distribution 
and various spatial distribution functions. The spatial 
distributions include the empty and occupied gap dis- 
tribution functions, and the density and velocity struc- 
ture functions. For the early time regime t ti, these 
quantities are numerically equivalent to the correspond- 
ing quantities of the sticky gas [2^ . In this paper, we 
focus on the late time regime t ^ ti, where there is de- 
viation from the sticky gas. 

Let Pi {p, t) be the probability that a box has mass 
density p at time t. In Fig. [21 the variation of Pi{p,t) 
with p is shown for different times. For small p {p < 10), 
Pi{p,t) increases with time to a nonzero constant. For 
intermediate p the slope increases to a constant. The 
cutoff increases to infinity with time. From these fea- 
tures, we conclude that when i ^ 1, Pi{p,t) approaches 
a non-zero time independent power law distribution, i.e. 

\\m lim Pi{p,t) ^ p~^K (3) 

Conservation of density ((p) = 1) implies that 71 > 2. 
Consistent with this, we find that the power law for 
the curves at various times extrapolate asymptotically 
to 7i ~ 2.83. The cutoff pmax{t) scales as C{t) ~ 
with (3 = 2/3. There are strong corrections to scaling. 
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FIG. 3: The variation of Pj*^'"' {x, t) with x is shown for times 
t = 16000, 32000, 64000 (the larger times have larger cutoffs). 
P2°\x,t) is shifted downwards for clarity. The straight lines 
have power —2.2. Inset: The plateaus of P2''\x,t) have the 
scaling form t''^^^ /{xt-^^^). 

as is clear from the change in apparent slope for different 
times. Hence, it is not possible to obtain data collapse by 
scaling unless one measures Pi {p, t) for even larger den- 
sities and times. The value of 71 depends on the value 
of a. For a varying from 3 to cxd, 71 varies from 2.30 to 
2.83. 

The data of Fig. [2] for our granular gas implies a 
scenario in which small density clusters do not get de- 
pleted from the system. This is in contrast to the 
sticky gas, where the probability of the mass clusters , 
Pis{m,t) ~ 'm~^/'^t~^/^^ for m <^ t}!'^ Q decays to zero 
with time. The coarsening process in the sticky gas is 
one of pure aggregation, transferring mass from smaller 
to larger mass scales. On the other hand, our model 
at late times effectively shows a combined aggregation 
and fragmentation dynamics. This simply comes as the 
"elastic" collisional break-ups at late times (and small 
velocities) compete with the aggregation due to inelas- 
ticity. The effective rates are such that mass loss and 
gain at small scales are balanced out. 

Further support to the above picture comes from the 
inter particle and inter hole gap distribution functions. 
Let P'^'°^ (x, t) be the probability of finding a gap of ex- 
actly X empty(e) or occupied(o) boxes. The variation of 
P2'''''\x,t) with X for different times is shown in Fig. [3l 
For small masses, they decay as a power law with power 
72 ~ —2.2. For large masses, P2°\x,t) has a plateau 
eventually cutoff at scales x ~ ^{t)- The shape of the 
plateau and cutoff is reminiscent of the the number dis- 
tribution of gaps in the sticky gas which have the form 
N2s{x,t) = t-^/^f2s{x/C{t)) with f2s{z) -> 1 for z < 1 
0. The large x of of Pg"-* (x, t) scales exactly as N2s{x, t) 
(see inset of Fig. Since / dxN2s ^ t~^^^, the area un- 
der the plateau will eventually go to zero. So while the 
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FIG. 4: The variation of Spp{k, t) with k is shown for different 
times. The inset shows the data scaled as in Eq. (|8]). The 
straight line has power —0.80. 



big gaps in our model scale as that of the sticky gas, the 
more important fact is that it has predominantly small 
gaps associated with the following power law: 

lim lim P^''-°\x,t) ^ x-"'^ (4) 

We have checked that the exponent 72 has no discernible 
dependence on a. 

The inter particle gap distribution has a bearing on the 
density-density and velocity-velocity correlation func- 
tions. In particular, the large power 72 suggests that 
the ordering process will be affected by the abundance of 
smaller gaps and Porod law could be violated. For the 
sticky gas, numerical results [2^, confirm that Porod 
law holds for both density-density and velocity-velocity 
correlation functions. We now discuss the case of our 
granular gas and show that it is different. 

Let Cpp(x,i) — {pi{t) pi+x[t)) be the equal time 
density-density correlation function. The structure func- 
tion Spp{k, t) is the Fourier transform of Cpp{x, t). Sim- 
ilarly we define the equal time velocity-velocity correla- 
tion function as Cvv{x,t) = {vi(t)vi^xit)), with its cor- 
responding structure function Smik,t). We show that 
the structure functions can be expressed in terms of the 
gap distribution through a mean field approximation. 
Cpp{x) is approximately equal to density square times 
the probability there is a non-zero density at x given 
there is a non-zero density at 0. Thus (setting p — 1), 
Cpp{x) « I]^o^'2n(a;), where P2n(x) is the probability 
of having exactly n empty gaps between and x with 
and X being occupied. Let Cpp{s) and P2''''^''(s) be the 
Laplace Transforms of Cpp (a;)and P2°''^\x) respectively. 
Approximating joint distributions by products of Individ- 
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FIG. 5: The variation of Swik, t) with k is shown for different 
times. The inset shows the data scaled as in Eq. (|9]). The 
straight line has power —0.80. 



ual distributions [231, we obtain 



P2n{s) 

Pais) 



1-p, 



(o) 



P. 



1 - P. 



(o) 



-, " > 1, (5) 
(6) 



where {x)o — J xP2°\x)dx. Eqs. ([5]) and (O give 



s {x),s^[l-Pi'\s)p!f\s)] 



(7) 



The correlation function Cyv{x) is exactly the same ex- 
cept for a factor of (w^). Thus, Cvv{s) ~ v^Cpp{s). 



Now, 1 - K 



(e,o) 



Also, 



£0 



since 



72 > 2. This implies that Cpp{s) ~ s^^~'^, where we have 
ignored the s^^ which contributes to 6{k) in the Fourier 
transform. Thus, Spp{k) will have the scaling form 



where f^ix) 



Spp{k)^C^-^-h{kC), (8) 

x'^" for a; 1, with 9p — 3 ~ j2- Using 
, we obtain 



where /4(x) ^ x~^^ for x 3> 1 with 9^ = 3 — 72. Since 
72 « 2.20, 9v — dp K, 0.80 which is very different from 
2.0 seen in Porod law in one dimension. 

The mean field approximation gives a good description 
of the actual problem. In Fig. [4] we show the variation of 
Spp{k, t) with k. The data collapses onto one curve when 
scaled as in Eq. ([8|) [see inset of Fig. [4j. The scahng 
function f^iz) varies as z""-^ for large z. In Fig. [H the 
variation of Svv{k,t) with k is shown. Again a good 
collapse is obtained when the data is scaled as in Eq. Q 



such that /4(z) ^ z "-^ for large z. We found that the 
exponents 9p and 9^ have no dependence on a. 

To summarize, we studied a model of freely cooling 
granular gas with velocity dependent restitution coeffi- 
cient. We showed the existence of a time scale ti ~ 
beyond which the system deviates from the sticky gas 
behavior. The effective dynamics in this regime is one of 
aggregation and fragmentation. As a result, the spatial 
distribution functions change their forms drastically. We 
found new power law exponents associated with the den- 
sity distribution function and empty and occupied gap 
distribution functions. The two-point spatial correlation 
functions violate Porod law and the structure function 
decay exponent is ~ 0.8 instead of the usual value 2. This 
deviation and the existence of power laws in the one-point 
functions indicate that the phase ordering is dominated 
by large scale fluctuations. We hope that experiments on 
cooling granular gases in quasi one-dimension would find 
the deviation in Porod law predicted by us. 
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